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We investigate the relevance of multiple-orbital and multipole effects during high-harmonic gen- 
eration (HHG). The time-dependent configuration-interaction singles (TDCIS) approach is used to 
study the impact of the detailed description of the residual electron-ion interaction on the HHG 
spectrum. We find that the shape and position of the Cooper minimum in the HHG spectrum of 
argon changes significantly whether or not interchannel interactions are taken into account. The 
HHG yield can be underestimated by up to 2 orders of magnitude in the energy regio of 30-50 eV. 
We show that the argument of low ionization probability is not sufficient to justify ignoring multiple- 
orbital contributions. Additionally, we find the HHG yield is sensitive to the nonspherical multipole 
character of the electron-ion interaction. 

PACS numbers: 32.80.RM,42.65.Ro,31.15.A- 



I. INTRODUCTION 

High-harmonic generation (HHG) is the key physical 
process underlying the generation of single attosecond 
pulses [1-3] and attosecond pulse trains [4-6], which are 
at the heart of attosecond science [7, 8]. In recent years, 
the rapid progress in HHG has led to applications rang- 
ing from atomic systems [9-11] over molecular systems 
[12, 13] to solid-state systems [14]. HHG has opened a 
new door to probe structural information [12, 15] as well 
electronic and nuclear dynamics [16-19] on fundamental 
time scales. 

The mechanism behind HHG is well captured in the 
three-step model [20], where in the first step the elec- 
tron is tunnel-ionized by a strong-field laser pulse, in the 
second step the electron is accelerated in the oscillating 
laser field, and finally in the last step the electron recol- 
lides with the parent ion and converts its excess energy 
into radiation energy in the extreme ultraviolet range 
[21, 22]. The maximum photon energy is given by the 
cut-ofi^ law, 1.32/p -|- 3.17C/p, where Ip is the ionization 
potential of the system and Up is the ponderomotive po- 
tential created by the intense laser field [20]. As a result, 
the heavier noble-gas atoms have lower cut-off energies 
than the lighter ones, whereas the HHG yield does in- 
crease with the atomic number [23, 24]. Previous works 
have shown that the recombination step can be directly 
related to photoionization [25, 26], enabling the retrieval 
of the electronic structure of the system [12, 25, 27, 28]. 
A strong focus has been, in particular, on molecular sys- 
tems [29-32]. 
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The two most common theoretical approaches for de- 
scribing HHG are the semiclassical strong-field approx- 
imation (SPA) [20], which has been extended to in- 
clude Coulomb-interaction corrections [25, 33-35], and 
the single-active-electron (SAE) approximation [36-39], 
where the electron-ion interaction for many-electron sys- 
tems is described by a model potential [37]. The SAE 
approach is computationally more demanding than the 
SPA approach, and, therefore, has been limited to atoms 
and systems like Hj. In the literature [25], the SAE 
approximation has often been referred to as solving the 
time-dependent Schrodinger equation (TDSE) . The SAE 
approach has some limitations. Por example, it ig- 
nores contributions from multiple orbitals. Intensive 
studies have been performed to understand the impact 
of multiple-orbital contributions in molecular systems, 
which are essential to understand the HHG spectrum and 
subsequently to extract electronic-structure information 
[40, 41]. Recently, it has been shown that even in atomic 
systems it is crucial to consider multiple-orbital effects 
[42]. 

In this paper, we investigate the importance of 
multiple-orbital (multichannel) contributions and mul- 
tipole effects in the residual electron-ion interaction on 
the Cooper minimum in the HHG spectrum of argon 
[37, 43, 44]. Both aspects are commonly ignored in SPA 
and SAE calculations. Multichannel interactions [45] go 
beyond the independent-particle picture and cannot be 
captured in the language of SPA and SAE, whereas mul- 
tipole effects could, in principle, arise even in a single- 
channel model such as SAE. However, in atomic SAE 
calculations, it is common to model the electron-ion in- 
teraction by a spherically symmetric potential [37]. The 
interaction of the liberated electron with the hole state 
(channel), from which it originates, is called intrachan- 
nel interaction and leads for large electron-ion distances 
to the 1/r behavior of the Coulomb potential. If the 
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liberated electron is influenced by other orbitals the in- 
teraction is called interchannel coupling [45] . The impor- 
tance of interchannel coupling for HHG has been shown 
for xenon, where a clear signature of the giant dipole 
resonance of the Ad subshell [45] known from photoion- 
ization studies has been directly observed in the HHG 
spectrum of xenon [42]. The theoretical model wc arc 
utilizing to capture these aspects is based on a time- 
dependent configuration-interaction singles (TDCIS) ap- 
proach [46, 47]. We have demonstrated in previous works 
that this TDCIS approach is ideal to study systemati- 
cally multichannel effects in situations involving ioniza- 
tion [48, 49]. 

The paper is organized as follows. In Sec. II, we give 
an overview of our theoretical method, which we use to 
systematically study the influence of various approxima- 
tions of the residual electron-ion interaction on the HHG 
spectrum. In Sec. Ill, we explain the system parameters 
used in our calculations. The results are discussed in de- 
tail in Sec. IV. Atomic units [50] are employed through- 
out, unless otherwise noted. 



II. THEORETICAL METHODS 

The time-dependent Schrodinger equation of an N- 
electron system exposed to a linearly polarized external 
electric field is given by 

^^^\^^<{t))^Hit)\^it)), (la) 
H{t) = Ho + Hi- E{t)z, (lb) 

where |5'(t)) is the full N-electron wave function and H{t) 
is the exact N-body Hamiltonian, which can be parti- 
tioned into three main parts: (1) Hq = F — irjW is 
the sum of the time-independent Fock operator F and a 
complex absorbing potential (CAP), where (x| W |x') = 
[r — rcAp]^ Q{r — rcAp) S{x — x') and 9(r) is the Hcavi- 
side step function. (2) The electron-electron interactions 
that cannot be described by the mean-field potential in 
Ho are captured by Hi (= Vc — Vhf — Ehf', for a de- 
tailed description of these quantities, see Ref. [47]). (3) 
The term E{t) z is the laser- matter interaction in the 
electric dipole approximation. The CAP serves a purely 
numerical purpose. It prevents artificial refiections of the 
ionized photoelectron from the radial grid boundary and 
is located far away from the atom. This is controlled by 
the parameter rcAP ■ 

Solving numerically the full N-electron system is cur- 
rently out of reach without making any approximations 
to the Hamiltonian or the wave function. In strong- 
field processes such as HHG, Eq. (la) is commonly re- 
duced to an effective one-electron system, where only 
one electron of the outermost valence shell is allowed 
to respond to the electric field and all other electrons 
are frozen or completely neglected. Here, we take an 
alternative way, by describing the full N-electron wave 



function and making no approximations to the Hamil- 
tonian. Specifically, we use the configuration-interaction 
language, where we assume the field-free ground state 
is well captured by the Hartree-Fock ground state ['I'o)- 
We consider only singly excited 1-particlc-l-hole config- 
urations (Iplh-configurations) The corresponding 
TDCIS N-clcctron wave function reads 

|*(t))=ao(t)|<I>o)+E""WI'5'")' (2a) 

i.a 

l*n = ^«At + <iea)l'i'o), (2b) 

where i, j and a, h refer to occupied orbitals and unoccu- 
pied (virtual) orbitals, respectively, in the Hartrec-Fock 
ground state ['I'o). The operators c)^ ^ and Ca.a create 
and annihilate, respectively, an electron in the orbital a 
with spin a. By restricting our wave function to Iplh- 
configurations, the interaction captured by Hi is the 
residual electron-ion interaction. The equations of mo- 
tion for the expansion coefficients ao(0 ^i^d a1{t) read: 

iao{t)^^E{t)Y,mim)o^1{t) (3a) 

■i,a 

zatit) = (6, - e,)a1{t) + ^ m\Hi I*?) «5(t) (3b) 

~E{t) ($,?|z|$o)ao(t) 
-i?(t)5^(4>»|z|$5)«5(t), 

b,3 

where tp are the orbital energies of the orbitals [y'p), 
which are eigenstates of the time-independet Fock oper- 
ator, i.e., Hq \lPp) = Cp \ipp). The expression (| stands for 
a dual vector with respect to the symmetric inner prod- 
uct, i.e., {(fp\(fq) = Sp,q, which differs from the Hcrmitian 
inner product. A detailed description of our implemen- 
tation of the TDCIS method can be found in Ref. [47]. 

The exact treatment of the residual electron-ion inter- 
action is numerically very demanding. In order to be able 
to treat the full electron-ion interaction, we are exploiting 
as much symmetry as possible. We have already used one 
symmetry with Eq. (2b), i.e., the total spin of the sys- 
tem (S ~ 0) is conserved. The second symmetry we are 
exploiting arises from the restriction to linearly polarized 
pulses and benefits us in two ways. Firstly, the orbital- 
angular-momentum projection TOq of the excited electron 
and the orbital-angular-momentum projection rrii of the 
hole state must be the same for each \^f). Secondly, 
the coefficients af{t) are the same whether an electron 
with orbital-angular-momcntum projection m or — m is 
excited. As a result, only the gerade parity configurations 
1$")^ need to be considered because ungerade parity con- 
figurations will not be populated due this symme- 
try. The gerade and ungerade parity configurations are 
defined as 

\'^^),/u--=^{\K^)±K'l))^ (4) 
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where the orbital indices ±a and ±i stand for triplets 
of quantum numbers (n, I, ±m) with n being the radial 
quantum number, I being the orbital angular momentum, 
and ±771 being the orbital-angular- momentum projection. 
The configuration |$^) with nii = = is a special 
case and has gerade parity. Since, as already mentioned. 



($^|z|$5) 



for linearly polarized light will not be populated 

and only 1$°)^ needs to be considered, we drop the index 
g such that \^f) refers to gerade parity configurations 
from now on. The matrix elements for the gerade parity 
configurations as they appear in Eqs. (3) are given by 



Z(+a,+b)^i.j - 2:(+j,+i)^a,b (5a) 

{y/2 , ttIq = = 
2 , ma = mj 7^ (5b) 
,ma^m^ 

(1 ,mj7^07^TOj, 

i ,m, = = mj, (5c) 

\/2 , otherwise, 



($o|i|<i>-) 
($-|Hi|<f>5) 



where we made use of the 



'^—a^—j\—i. 

lations m 



symmetries .^j) 

-b,-i and V^a,+j;+i,+b 



i)i'>^+a,+j;+b.+i — I'-a, 

= w+a.-j;+i,-6- In the last equation the re- 
= m.i and m?, = nij are used, which holds only 
in the case of linearly polarized light. The round paren- 
theses in the indices of the matrix elements indicate the 
symmetric inner product mentioned above [47]. 

In the following, we study in detail three scenarios 
for Hi: (1) no approximation is made and the residual 
Coulomb interaction is treated exactly within the CIS 
configuration space, (2) only intrachannel interactions 
[{^f\Hi = if i ^ j] are considered, (3) a sym- 
metrized version of the intrachannel interaction is used 
such that the angular momentum of the excited electron 
cannot be changed, thus simulating a spherically sym- 
metric ion potential. When only intrachannel interac- 
tions arc allowed, different orbitals will behave almost in- 
dependently. Only via the ground-state depopulation can 
they indirectly influence each other. The symmetrization 
in model (3) is done by averaging over all hole states 
within each (n, I) subshell such that the excited electron 
sees only a spherically symmetric ion. The symmetrized 
matrix elements read 



symm 



(a,i;i,b) ■ 21 



'^-p^'^V{a,t;b,i)^ (6a) 
^-^'^V{a,^■^,b)^ (6b) 



where in Eq. (6b) we additionally set rua = rrii — rrib 
before we perform the sum. This step can be justified, 
since we are using linearly polarized light and our model 
can only have 1-particle-l-hole configurations with TOq = 
rrii. In both cases, one finds that the symmetrized matrix 
elements are proportional to (5;^,;^ and i5m„,m6- 

The HHG spectrum, which is calculated via the expec- 
tation value of the electric dipole moment (z) (t), reads 



[51, 52] 

S{lo) 
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J — OO 
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Next to the HHG spectrum, we will focus our discussion 
also on the hole populations Pi{t) generated during the 
HHG process. These populations are calculated with the 
help of the ion density matrix, which is described in detail 
in Ref. [47]. The ground-state population is given by 
Po(t) = l«o(i)l'- 



III. NUMERICAL DISCUSSION 

All presented results were calculated with the XCID 
package [53], which makes explicit use of the symme- 
tries discussed in Sec. II. All argon calculations pre- 
sented in Sec. IV were done for a laser pulse with a 
peak field strength of E^max = 0.085, a carrier frequency 
uj = 0.057 (« 800 nm), and a FWHM pulse duration of 
T = 413 (w 10 fs). The classical turning radius of the 
electron for such a pulse is thhc = -E'max/'^^ ~ 26. As 
described in Ref. [47], we use a non- uniform grid with 
the mapping function 



■2 1-x + C 



,xe[-i,i]. 



(8) 



All calculations were done with a radial grid radius 
''max ~ 120, 480 radial grid points, and mapping parame- 
ter C = 1.0. The CAP starts at a radius rcAP = 100 and 
has a strength rj = 0.01. The maximum angular momen- 
tum employed was /max = 80. Furthermore, we find that 
for excited electrons with an orbital angular momentum 
I > 6 the multipole terms with Lc > are negligibly 
small and the dominant Hi contribution comes from the 
monopole term. It is, therefore, a good approximation to 
consider only the monopole term of Hi when any orbital 
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angular momentum of the involved orbtials is larger than 
6 [for details see Ref. [47]]. 



IV. RESULTS 

We begin our discussion with the single-channel model 
by allowing only the 3po orbital to be active. In Fig. 
1, we compare the HHG spectra and the depopula- 
tions of the ground state of argon for a spherically sym- 
metric electron-ion interaction (labeled symmetric) with 
the exact electron- ion interaction (labeled intrachannel) . 
Note, interchannel contributions do not exist in a single- 
channel model. The spherically symmetric Hi has no 
tensorial multipole moments besides a monopole term, 
since the angular momentum of the electron cannot be 
changed. In the intrachannel and interchannel models all 
multipole contributions are included in Hi. The depop- 
ulations [shown in Fig. 1 (b)] show only small deviations 
during and after the pulse. The final depopulation prob- 
abilities are almost identical. Similarly, the HHG spectra 
[see Fig. 1 (a)] show only small differences in the energy 
region of 30-50 cV, where also the photoionization cross 
sections (not shown) differ by up to 30% from each other. 
The Cooper minimum in the HHG spectra can be repro- 
duced and lies between 40-50 eV. The low curvature of 
the shape of the Cooper minimum prevents a more pre- 
cise localization of the minimum. 

Due to the costly treatment of the residual electron-ion 
interaction it is common to reduce an HHG calculation 
to a one-electron calculation, where the electron moves in 
a local, spherically symmetric model potential, which de- 
scribes the correct behavior for short and long distances 
of the electron-ion interaction and reproduces the ion- 
ization potential. This is the SAE approach, which is 
in spirit very close to our single-channel model with a 
spherically symmetric Hi (green dashed line in Fig. 1). 
However, there exists one major difference to typical SAE 
calculations. In the SAE approach one only describes one 
electron, which can move freely everywhere on the pscu- 
dopotential surface, and does not fulfill the Pauli princi- 
ple, meaning there is no mechanism in this approach that 
can prevent the electron to move into orbitals that are 
already occupied by the A^ — 1 frozen electrons. Enforcing 
the Pauli principle is critical for the one-electron reduced 
density matrix to be N-representable, that is, to represent 
a realistic N-clcctron system [54, 55]. In our theory, we 
describe always the entire A'-electron wave function and 
due to the anticommutator relation of the creation and 
annihilation operators in Eq. (2b) the Pauli principle is 
ensured at all times. Recent works for molecular systems 
have pointed out the importance of the Pauli principle 
particularly for tomographic purposes [56, 57]. 

Now we consider the impact of the approximation of 
the residual electron-ion interaction in the multichannel 
scenario, where we allow also <ip±i electrons and the 3s 
electrons to get ionized. All differences seen in Figs. 2, 
3, and 4 originate solely from physics within the 3p man- 
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FIG. 1. (color online) The HHG spectra (a) and the depopu- 
lations of the ground state (b) of argon are compared for the 
intrachannel approximation (red sohd line) and for the sym- 
metrized intrachannel approximation (green dashed line). In 
both cases only 3po (single-channel) is active. The intensity 
profile of the pulse is shown in (b). The pulse parameters are: 
Smax = 0.085, cj = 0.057 (^ 800 nm), and r = 413 (^ 10 fs). 



ifold. The impact of the 3s orbital is rather small on 
the HHG spectra as well as on the depopulation of the 
ground state. That is not a surprise due to the high 3s 
ionization potential, which is around 18 cV higher than 
the ionization potential of the 3p orbitals. The hole pop- 
ulation of 3s (not shown) is over 100 times smaller than 
the hole populations in the 3p shell. 

In Fig. 2 the HHG spectra and depopulations are 
shown for different approximations of Hi. The simpli- 
fication to a spherically symmetric potential [see Fig. 2 
(a,b)] underestimates the HHG spectrum by up to two 
order of magnitude in the energy region of 30-50 eV. In 
addition, the shape of the Cooper minimum has now 
drastically changed. The position of the Cooper mini- 
mum is much more clearly defined in the interchannel 
case and lies slightly above 50 eV as found in recent ex- 
periments [37, 43, 44]. For the symmetrized intrachannel 
approximation, the Cooper minimum lies between 40- 
50 eV, similar to the single- channel results. The depop- 
ulation dynamics is not affected by the approximation of 
the electron-ion interaction. In Fig. 2(c,d) the results 
of the intrachannel calculation are compared with the 
results obtained from the symmetrized intrachannel ap- 
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FIG. 2. (color online) The HHG spectra (a,c) and the depopulations of the ground state (b,d) of argon are shown for 3s and 
all 3p orbitals active. In (a,b) the interchannel and in (c,d) the intrachannel approximation is compared with the symmetrized 
intrachannel approximation. The pulse parameters are the same as in Fig. 1. 



proximation. The relative differences between these two 
models never exceed a factor larger than 2 and are con- 
fined to the energy region of 30-50 eV. The origin of these 
differences is the lack of multipole effects in the sym- 
metrized intrachannel model. In comparison to the in- 
terchannel results, the multipole effects are much smaller 
than the interchannel effects seen in Fig. 2(a). 

The photon energy range of 30-50 eV corresponds to 
a rccollision electron energy range of 15-35 eV and a de 
Broglie wavelength of 2-3 A. It seems that electrons with 
these wavelengths arc most sensitive to the exact residual 
ion-electron interaction and, therefore, simplifications of 
the interaction become most evident in the correspond- 
ing photon energy regime. In the same energy region 
the photoionization cross section is most sensitive to the 
approximation made to the electron- ion interaction [45]. 



Our calculations (not shown) confirm that the differences 
in the cross sections between the intrachannel and sym- 
metrized intrachannel model are quite small, where as the 
differences to the interchannel model are up to one order 
of magnitude larger (and can reach values up to 20 Mb) . 
The fact that the photoionization cross sections and the 
HHG spectra behave similarly (for different approxima- 
tions to the electron-ion interaction) supports the picture 
that HHG has a close connection to photoionization. In 
the limit that the electron de Broglie wavelength is much 
longer or shorter than the characteristic length scale of 
the residual electron-ion interaction, the results do not 
depend on the detailed structure of electron-ion interac- 
tion. This may explain why the HHG spectrum does not 
alter significantly for photon energies smaller than 20 eV 
and photon energies close to the cut-off region. 
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FIG. 3. (color online) The HHG spectra (a.c) and the depopulations of the ground state (b,d) of argon are shown. The single- 
channel (only 3po active) calculation is compared with the multichannel calculation (3s and all 3p active). No approximation 
(a,b) and the symmetrized intrachannel approximation (c,d) are made to the electron-ion interaction. The pulse parameters 
are the same as in Fig. 1. 



We have seen that approximations to the electron-ion 
interaction can cause differences in the HHG yield, par- 
ticularly, when multiple orbitals are considered. Not 
all electron-ion approximations make the HHG spectrum 
sensitive to whether a single orbital or multiple orbitals 
are active. In Fig. 3 we compare the HHG spectrum 
and the depopulation of the ground state for a single- 
channel (only 3po active) and for a multichannel (all 3p 
and 3s arc active) calculations. Figures 3(a-b) are cal- 
culated with the exact Hi term including interchannel 
and multipole effects in the residual electron-ion inter- 
action. The HHG signal strength for the single-channel 
calculation is strongly reduced in the spectral range of 
30-50 eV, whereas the signal is slightly enhanced in the 
cut-off region. When the symmetrized intrachannel ap- 



proximation is made [sec Figs. 3(c-d)], the HHG spectra 
are almost identical whether single-channel or multichan- 
nel calculations are performed. This stands in contrast 
to the interchannel results, where the interchannel cou- 
pling causes strong differences in the HHG yield. The 
depopulation is overestimated by sa 25% regardless of 
the approximation made to the electron-ion interaction. 

We have seen the strong differences in the HHG spec- 
tra and in the depopulations depending on whether only 
3po or all 3p electrons are active (the contributions from 
3s are negligible small). Does that also mean the ion- 
ization probabilities of 3p±i are comparable with 3po ? 
The ratio {psp^ + P3p_i)/p3po is shown in Fig. 4 for the 
differenet models. Note that for linearly polarized light 
P3p+i — P3p_i- Before the pulse, all ratios are close to 
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FIG. 4. (color online) The ratio of the hole populations p-^p^ + 
P3p_i and p3pg is shown for different approximations to the 
electron-ion interaction. The pulse parameters are the same 
as in Fig. 1. 



0.8 and drop to « 0.1 after the pulse. Both intrachannel 
models lead even to the same final ratio. The small ra- 
tios after the pulse show that, at least ultimately, mainly 
the 3po orbital gets ionized, which is at the heart of the 
SAE approximation. The oscillations in the ratios during 
the pulse are in phase with the oscillations in the electric 
field. They are a direct consequence of the projection of 
the wave function onto the field-free states in the pres- 
ence of the laser field. In contrast to the HHG spectra 
(see Fig. 2), the effects of the intrachannel or symmet- 
ric approximation on the population dynamics arc quite 
small. 



V. CONCLUSION 

We have described the HHG process with a many-body 
approach, namely TDCIS, where we describe the entire 
N-electron wave function. This allows us to fulfill the 
Pauli principle at all times. Our results show that mul- 
tichannel effects in the residual electron-ion interaction, 
which is a combination of the bare nuclear potential and 
the electron-electron interaction for many-electron sys- 
tems, have a significant influence on the HHG spectrum. 
They cannot generally be neglected for atoms and specifi- 
cally not for molecules as recent experiments have shown 
[40, 41]. We have demonstrated that orbitals, despite 
their relatively low ionization probability by the end of 



the pulse, can lead to surprisingly large modifications 
of up to 2 orders of magnitude in the HHG spectrum 
(especially in the energy region of 30-50 cV). While we 
confirm that after the end of the pulse, the populations 
of the 3p±i orbitals are relatively small, their contribu- 
tions during the pulse are not small and have indirectly 
through interchannel coupling a significant impact on the 
HHG yield. 

We saw that neglegting interchannel interactions lead 
to large changes in the HHG yield. Multipole effects in- 
fluence the spectra but not as dramatically as interchan- 
nel effects do. All deviations in the HHG yield are in the 
30-50 cV energy region, which corresponds to a dc Broglie 
wavelength of the recolliding electron between 2-3 A. This 
coincides with the characteristic length scale on which the 
electron-ion interaction goes over into a pure long-range 
1/r potential. In contrast to the large disagreement in 
the HHG spectra between the single-channel and multi- 
channel calculations including interchannel interactions 
we found that by using the symmetric interaction the 
HHG spectra look quite the same whether or not a sin- 
gle or multiple channels participate in the HHG process. 
This comparison directly shows that the population of 
an orbital does not map directly to its importance in the 
HHG mechanism. 

All these observations demonstrate that many-body 
effects enter in the HHG spectrum and need to be un- 
derstood in order to successfully use them for tomo- 
graphic imaging [12]. The time-dependent configuration- 
interaction approach provides a clear pathway how these 
and higher-order effects can be taken into account. Re- 
cent works [58] have suggested that multiclcctron excita- 
tions are not a dominant factor. All essential multielec- 
tron effects can be captured by single-electron excitations 
including interchannel interactions. This makes the TD- 
CIS approach perfectly suited for studying many-body 
effects in HHG. 
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